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A new method is discussed that may be used to correct wind tunnel strain 
gage balance load predictions for the influence of residual temperature effects 
at the location of the strain gages. The method was designed for the iterative 
analysis technique that is used in the aerospace testing community to predict bal- 
ance loads from strain— gage outputs during a wind tunnel test. The new method 
implicitly applies temperature corrections to the gage outputs during the load 
iteration process. Therefore, it can use uncorrected gage outputs directly as 
input for the load calculations. The new method is applied in several steps. 
First, balance calibration data is analyzed in the usual manner assuming that 
the balance temperature was kept constant during the calibration. Then, the 
temperature difference relative to the calibration temperature is introduced as a 
new independent variable for each strain— gage output. Therefore, sensors must 
exist near the strain— gages so that the required temperature differences can be 
measured during the wind tunnel test. In addition, the format of the regression 
coefficient matrix needs to be extended so that it can support the new inde- 
pendent variables. In the next step, the extended regression coefficient matrix 
of the original calibration data is modified by using the manufacturer specified 
temperature sensitivity of each strain— gage as the regression coefficient of the 
corresponding temperature difference variable. Finally, the modified regression 
coefficient matrix is converted to a data reduction matrix that the iterative 
analysis technique needs for the calculation of balance loads. Original calibra- 
tion data and modified check load data of NASA’s MC60D balance are used to 
illustrate the new method. 
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= axial force component 

= parts of data reduction matrix that are used with the primary load iteration equation 
= balance load 

= vector of predicted balance loads for load iteration step £ 

= load dependent matrix that is used in the load iteration equation 
m strain-gage index 
= total number of strain-gage outputs 
= total number of balance load components 
= forward normal force component 
= aft normal force component 
= electrical output of a strain-gage 
= electrical outputs of a strain-gages 
= rolling moment 
= forward side force component 
= aft side force component 
= temperature at strain-gage location 

= temperature difference at strain-gage location as dependent variable 
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Tli, T2i, ... = temperature difference at strain-gage location as independent variable 

T 0 = uniform calibration temperature of a strain-gage balance 


ao, an, ... 

AR 

AT 

£ 


= regression model coefficients of fitted gage outputs 

= vector of strain-gage output differences that is used for iterative load calculation 
= difference between temperature at strain-gage and calibration temperature of balance 
= load iteration step index 


I. Introduction 

Different techniques are used in the aerospace testing community to predict strain-gage balance loads 
from measured gage outputs during a wind tunnel test. The load prediction is usually based on the result of 
a multivariate regression analysis of strain-gage balance calibration data. The iterative analysis technique , 
for example, first fits strain-gage outputs of a balance as a function of the loads that were applied during 
the calibration. Then, the regression coefficients are converted to data reduction matrix coefficients so that 
loads can be predicted from measured gage outputs using a load iteration scheme (see Ref. [1] for a detailed 
description of the iterative analysis technique ). 

By design, the iterative analysis technique makes the assumption that the number of independent vari- 
ables (i.e., balance loads) must equal the number of dependent variables (i.e., strain-gage outputs). This 
limitation made it difficult in the past to include corrections for temperature effects directly in the regression 
model of a strain-gage output. 

In principle, corrections for temperature effects may be estimated by using the temperature sensitivity 
of a strain-gage in combination with a new independent variable that is equal to the difference between 
the temperature of the balance during a wind tunnel test and the calibration temperature of the balance. 
Recently, an extension of the iterative analysis technique was introduced that allows an analyst to use more 
independent than dependent variables for the analysis of balance calibration data (for a description of the 
extension of the iterative analysis technique see Ref. [2]). This improvement makes it now possible to directly 
include temperature corrections during the load iteration process. 

In general, temperature effects on the gage outputs may be included in the balance load calculation 
process using one of two processing approaches. The first approach corrects the gage outputs for temperature 
effects before the load iteration takes place. This approach has traditionally been used in the aerospace testing 
community. The second approach , i.e., the new approach, applies the temperature correction during the load 
iteration process by extending and editing the original regression model of the gage outputs. This leads to a 
new data reduction matrix that may be used with the load iteration scheme of the iterative analysis technique 
to predict the balance loads. 

In the next section of the paper the two approaches are explained in more detail. First, the theoretical 
background of the two approaches is reviewed. Then, a realistic simulation of a temperature sensitive balance 
check load data set is discussed to illustrate basic elements of the newly introduced second approach. 


II. Inclusion of Temperature Sensitivity in Iterative Analysis Technique 

A. General Remarks 

Sometimes the accuracy of iterative strain-gage balance load predictions may be improved by including a 
residual temperature correction in the regression model of the gage outputs. This residual correction ignores 
temperature gradients that may exist inside of the balance. Therefore, it only corrects for the difference 
between the uniform balance calibration temperature and the temperature of the balance during the wind 
tunnel test. 

The temperature correction of the gage outputs may be a simple linear correction that uses (i) the 
temperature sensitivity and (ii) the temperature difference at the strain-gage as input. The first input, 
i.e., the temperature sensitivity of a strain-gage, is defined as the derivative of the electrical output of a 
strain-gage with respect to the temperature at or near the gage location. Then, we can write 


temperature sensitivity 


d R 

Iff 


( k ) 
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where R is the electrical output of the gage, T is the temperature at or near the gage, and k is the gage 
index. A numerical value of the derivative may be obtained by performing a sensitivity experiment for each 
strain-gage. Alternatively, the derivative may be supplied by the manufacturer of the strain-gage. 

Now, a simple way needs to be found so that the temperature sensitivity can be included in the regression 
model of the gage outputs that the iterative analysis technique uses. Therefore, the regression model of the 
gage outputs at a constant calibration temperature T C) needs to be revisited. It is given by the following 
equation (for more detail see Ref. [1], Eq. (3.1.3), or, Ref. [3], Eq. (13)): 


R(k,T 0 ) — cko(/c,T 0 ) + aq(fc, T 0 ) • T(l) T (^(/c, T 0 ) • _F(2) T ••• 

intercept linear and non-linear terms 


1 < k < l (1) 


The original set of independent variables and regression coefficients used in Eq. (1) can be summarized 
as follows: 


original set of independent variables ==> F( 1), F( 2), •••, F(n ) 
original set of regression coefficients ao (fc, T 0 ), ai(fc, T 0 ), ••• 


A linear temperature correction of the gage outputs may be included in Eq. (1) by simply adding the 
product of the temperature sensitivity and the temperature difference relative to the calibration temperature 
to the left hand side of Eq. (1). Then, we get: 


where 


R(k,T) = R(k,T 0 ) + 


d r 
d t 


(k) 


A T{k) 


1 < k < l 


temp, sensitivity 


A T(k) = T(k) 


(2a) 


( 26 ) 


At this point a user of the iterative analysis technique may use one of two processing approaches in 
order to include temperature corrections. The first approach explicitly applies the temperature correction to 
the gage outputs before the load iteration process. The second approach implicitly applies the temperature 
correction to the gage outputs during the load iteration process. The two processing approaches are explained 
in more detail in the following sections. 


B. Explicit Application of Temperature Corrections 

This first processing approach is traditionally applied in the aerospace testing community. It corrects 
the measured gage outputs for the difference between the calibration temperature and the temperature of 
the balance during the wind tunnel test before the load iteration begins. The required relationship between 
the gage output at the calibration temperature T 0 and the measured gage output at the balance temperature 
T during the wind tunnel test is obtained after rearranging Eq. (2a). Then, we get 


R(k,T 0 


corr. value 


R(k, T) 


meas. value 


d R 
d T 


(k) 


A T{k) 


1 < k < l 


temp. corr. 


( 3 ) 


where the temperature difference AT is given by Eq. (2b). It is important to point out that the corrected 
gage output R(k, T 0 ) and not the measured gage output R(k, T ) must be used as input for the balance load 
calculation because the original data reduction matrix is only valid for the calibration temperature T 0 . The 
application of Eq. (3) also has the advantage that no change of the original data reduction matrix is required. 
The gage outputs are simply transformed back to the electrical output T(/c, T 0 ) that the balance would have 
had if the temperature during the wind tunnel test would have matched the original calibration temperature. 


C. Implicit Application of Temperature Corrections 

This second processing approach implicitly applies the temperature correction to the gage outputs during 
the load iteration process. It has the advantage that the temperature effects are not treated separately from 
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the load calculation. Instead, they are directly included in the load iteration process. The second approach, 
however, has the disadvantage that the original regression analysis result needs to be modified so that the 
temperature difference A T(k) (see Eq. (2b)) may be used as an independent variable for the load iteration 
process. This modification may be better understood by reviewing the regression model of the gage outputs 
that the iterative analysis technique uses. After using the right hand side of Eq. (1) to replace the gage 
output R(k,T c> ) in Eq. (2a), we get the following expression: 


R(k, T) — ao (fc, T 0 ) + aq(/c, T 0 ) • F(l) + ct 2 (k, T 0 ) • F(2) + ••• + 

' V ' 

regression model of data recorded at calibration temperature T 0 

Equation (4a) can also be interpreted in a different way. We can write: 

ind. var. ind. var. 

R(k, T) = ao(fc, T 0 ) + oq (k,T 0 ) • F(l) + ol^ (k,T 0 ) • F( 2) + ••• + 

coef. coef. coef. 


■§£<*> 


• A T{k) (4a) 


linear correction 


ind. var. 


d R 
d T 


(k) 


• A T(k) (46) 


coef. 


Consequently, the extended set of independent variables and regression coefficients used in Eq. (46) can 
be summarized as follows: 


extended set of independent variables ==> F( 1), F( 2), •••, F(n), A T(k) 

d R 

extended set of regression coefficients ==> ao(fc, T ), ai(fc, T 0 ), •••, (fc) 

The regresssion coefficients ao(/c, T 0 ), aq(/c, T 0 ), ... of Eq. (46) are the result of the regression analysis 
of the original balance calibration data that was recorded at calibration temperature T 0 . The additional 
regression coefficient dR/dT is the manufacturer specified temperature sensitivity of the gage. 

It remains to include the additional independent variable A T(k) in the load iteration process that the 
iterative analysis technique uses. This additional independent variable needs to be added to the original 
data reduction matrix using an approach that was first introduced in Ref. [2] . At this point only basic steps 
of the preparation of the new data reduction matrix are summarized. They will be explained in more detail 
using the data example that is described in the next section. 

The flowchart shown in Fig. 1 summarizes the most important steps that are needed for the prepa- 
ration of the extended data reduction matrix that takes residual temperature effects into account. First, 
the balance is calibrated in the usual manner at a constant calibration temperature T Q ( Step 1 in Fig. 1). 
Then, the global regression solution of the original balance calibration data is found. This initial regression 
analysis solution is given by the regression coefficients ao(k. T ), oq (k:. T ), ... that are used in Eqs. (1) and 
(46) above. The coefficients of all strain-gages of the balance are assembled in the regression coefficient 
matrix after the original regression analysis has been completed ( Step 2 in Fig. 1). Now, the temperature 
sensitivities of the balance gages are determined in a separate sensitivity experiment ( Step 3 in Fig. 1). 
Alternatively, the sensitivities may have been supplied by the manufacturer of the strain-gages. In the next 
step, the original regression coefficient matrix is modified ( Step 4 in Fig. 1). Therefore, the format of the 
regression coefficient matrix needs to be extended so that it supports the additional independent variable set. 
For a six-component balance, for example, this means that the regression coefficient matrix is transformed 
from the original 6 x 97 format to the new 12 x 337 format because each strain-gage output is accompa- 
nied by a corresponding temperature difference measurement. Then, the temperature sensitivity of each 
strain-gage is entered as the linear coefficient of the corresponding temperature difference variable (i.e., of 
the new independent variable ). In addition, using ideas presented in Ref. [2], a value of “1.0” is used as the 
coefficient of the new independent variable in the regression model of the matching new dependent variable . 
Finally, the modified regression coefficient matrix in the new 12 x 337 format is transformed to the corre- 
sponding data reduction matrix that contains matrices needed for the load iteration process. 
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Data from the calibration of a six-component force balance and a realistic simulation of check load data 
are used in the next section of the paper to illustrate the five steps that were discussed above. 

III. Discussion of Example 

Manual calibration data of NASA’s MC60D balance was selected to illustrate the newly introduced 
implicit processing approach that was discussed in the previous section. The MC60 balance was manufactured 
and calibrated by Triumph/Force Measurement Systems in 2008. Table 1 below summarizes important 
features of the balance and the calibration data set that was used for the present study. 

Table 1: Balance and calibration data set characteristics of the MC60D balance. 


CHARACTERISTIC 

DESCRIPTION 

BALANCE TYPE 

force balance 

DIAMETER 

2.0 [in] 

EXCITATION VOLTAGE 

6.0 [V] 

CALIBRATION DATE 

December 2008 

CALIBRATION METHOD 

manual calibration 

NUMBER OF CALIBRATION POINTS 

175 


Table 2 lists load capacities of the balance: 

Table 2: Load capacities of the MC60D balance. 



Nl, lbs 

N2, lbs 

SI, lbs 

S2, lbs 

RM, in-lbs 

AF, lbs 

CAPACITY 

2500 

2500 

1250 

1250 

5000 

700 


The balance was calibrated using a manual calibration approach. A total of 175 calibration points were 
recorded in 17 load series. The temperature was kept constant during the entire calibration. For the present 
study it was assumed that the calibration data was recorded at a nominal temperature of 70 [degF]. 

In the next step, the original balance calibration data was analyzed using the iterative analysis technique . 
Afterwards, the regression coefficient matrix in extended format (six gage outputs plus six temperature 
differences) was sent to a text file for further processing. Figure 2a shows the regression coefficient matrix 
in the extended format before its modification. The matrix is in the extended 12 x 337 format (only 50 of 
337 rows are shown). The green color highlights the regression coefficients of the first gage output (Rl) that 
are the result of the analysis of the original calibration data. 

Figure 2b shows the regression coefficient matrix after it was modified. The modified regression coeffi- 
cients are highlighted in the blue box. Two modifications may be seen by comparing Fig. 2a with Fig. 2b. 
Modification 1: The red numbers are simulated temperature sensitivities of the six balance gages that were 
entered in the original matrix. These simulated sensitivities are much larger than the true temperature 
sensitivities of the six gages of the MC60D balance. They were only artificially increased for the purpose of 
this study to better illustrate the application of the new method. Modification 2: A value of “1.0” has to be 
used as the coefficient of a new dependent variable that matches a corresponding independent variable (see, 
e.g., the coefficient of the independent variable Tld). 

Finally, the modified regression coefficient matrix (Fig. 2b) was converted to the corresponding extended 
data reduction matrix. It is assumed that this data reduction matrix supports the following load iteration 
equation (from Ref. [3], p. 12, Table 4): 


F 



AR - 


Ci'Ca 

matrix 


• H(F C _0 


( 5 ) 
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Equation (5) is the load iteration equation that is recommended in Ref. [1]. It is also derived in great 
detail in Ref. [3]. The equation parts C x and C x C 2 are contained in the data reduction matrix (see 
also Fig. 3). The conversion from regression coefficient matrix to data reduction matrix is a well defined 
mathematical operation that is described in Ref. [1] and Ref. [3]. Figure 3 shows the converted data reduction 
matrix that was obtained for the modified regression coefficient matrix that is shown in Fig. 2b. The blue 
data block shows terms of the data reduction matrix that are related to the temperature sensitivities and 
the temperature differences Tli, ..., T6i that were introduced as new set of independent variables. 

The modified data reduction matrix shown in Fig. 3 was tested using a set of simulated check loads 
of the MC60D balance. The original check loads remained unchanged. The gage outputs of the check 
loads, however, were modified by adding a gage output change to the original gage output. The output 
change was computed by using the simulated temperature sensitivities depicted in Fig. 2b and an assumed 
temperature difference between the actual and calibration temperature at the location of the strain-gage as 
input. The temperature differences were randomly assigned using values between -20 degF and +40 degF. 
First, the check loads were processed using the gage outputs of the check load points and assuming that 
all temperature differences of the check loads were zero. The corresponding check load residuals are shown 
in Fig. 4a. Then, the check loads were processed using both the gage outputs and the actual values of the 
temperature differences as input. The corresponding check load residuals are shown in Fig. 4b. Comparing 
the check load residuals of the six balance loads components without and with temperature correction we see 
that the residuals are noticable smaller when the temperature differences were used. This result illustrates 
that the implicit use of temperature corrections in an extended data reduction matrix works as intended. 

IV. Summary and Conclusions 

A new method was discussed that may be used to include residual temperature corrections of 
strain-gage outputs in balance load calculations. The new method was specifically developed for the 
iterative analysis technique . It has the advantage that it leads to a data reduction matrix that is also a 
function of the temperature of the balance during a wind tunnel test. Therefore, the application of residual 
temperature corrections no longer needs to be separated from the load iteration process. 

The new method implicitly applies the temperature correction during the load iteration by processing 
an extended data reduction matrix that uses the temperature difference between the balance temperature 
during the wind tunnel test and the balance temperature during the calibration as a new input variable. 
Manufacturer specified temperature sensitivities of strain-gages are directly used as regression coefficients 
of the related temperature difference variables. Calibration data and simulated check load data of NASA’s 
MC60D six-component strain-gage balance were used to successfully test the proposed new method. The 
analysis of check load residuals, i.e., of the difference between applied and predicted check loads, showed 
that the direct inclusion of temperature sensitivities of the strain-gages in the data reduction matrix kept 
temperature dependent load prediction errors to a minimum. 

The new method applies residual temperature corrections during the load iteration process as the 
temperature sensitivies are hidden in the data reduction matrix coefficients. Therefore, gage outputs may 
directly be used as input for the load calculation that are not corrected for residual temperature effects. It 
is important, however, to point out that the new approach takes advantage of unique capabilities of NASA’s 
BALFIT software (see Ref. [4] for a basic description of BALFIT). Since 2010 BALFIT can process a data 
reduction matrix with an extended set of independent variables (see Ref. [2]). In addition, the inclusion of 
temperature sensitivities in the data reduction matrix has recently been automated. Therefore, as far as 
BALFIT ’s user is concerned, the effort required to prepare a data reduction without temperature correction 
is essentially identical to the effort required to prepare a data reduction matrix with temperature correction. 
Temperature sensitivities only have to be specified in BALFIT ’s calibration data input file whenever a 
temperature dependent data reduction matrix has to be generated. Then, BALFIT executes steps 2, 4, and 
5 shown in Fig. 1 without requiring any user interaction. Most existing balance calibration data analysis 
software packages do not have these capabilities. In those cases, the residual temperature correction has to 
be applied explicitly to the gage outputs before the load iteration starts (see Eq. (3) ) as the unmodified data 
reduction matrix is only valid for the original calibration temperature of the balance. 
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STEP 1 BALANCE CALIBRATION AT CONSTANT TEMPERATURE 
CALIBRATE BALANCE AS USUAL AT A UNIFORM TEMPERATURE (T 0 ) 
(BALANCE SHOULD NOT EXPERIENCE A TEMPERATURE GRADIENT) 


S TEP 3 SENSITIVITY EXPERIMENT 
CONDUCT EXPERIMENT TO GET 
TEMPERATURE SENSITIVITY 
OF EACH STRAIN-GAGE 


APPLIED LOADS, F(j) 
MEASURED GAGE OUTPUTS, R(k) 
CALIBRATION TEMPERATURE, T 0 


J 


TEMPERATURE 
SENSITIVITIES 
dR(k) /dT(k) 


STEP 2 REGRESSION COEFFICIENT MATRIX CALCULATION 
FIT GAGE OUTPUTS AS A FUNCTION OF LOADS 
AND GET REGRESSION COEFFICIENT MATRIX 




i 




STEP 4 REGRESSION COEFFICIENT MATRIX MODIFICATION 
INTRODUCE TEMPERATURE DIFFERENCE (AT) AT EACH GAGE AS A NEW 
INDEPENDENT VARIABLE AND USE TEMPERATURE SENSITIVITIES AS 
REGRESSION COEFFICIENTS OF THE TEMPERATURE DIFFERENCE TERMS 
IN A MODIFIED VERSION OF THE REGRESSION COEFFICIENT MATRIX 


STEP 5 DATA REDUCTION MATRIX CALCULATION 
CONVERT MODIFIED REGRESSION COEFFICIENT MATRIX TO THE FINAL 
DATA REDUCTION MATRIX ; THIS MATRIX WILL USE (I) BALANCE GAGE 
OUTPUTS AND (II) TEMPERATURE DIFFERENCES AT EACH GAGE AS INPUT 
FOR THE CALCULATION OF BALANCE LOADS DURING A WIND TUNNEL TEST 



Fig. 1 Preparation of a data reduction matrix that uses temperature sensitivities. 
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Fig. 2a Regression coefficient matrix before modification (50 of 337 rows shown). 
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Fig. 2b Regression coefficient matrix after modification (50 of 337 rows shown). 
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Fig. 3 Data reduction matrix of modified regression coefficient matrix (37 of 349 rows shown). 
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ARITH. MEAN = -1.0876 [lbs] (-0.0870 %) ; STD. DEV. = 2.9507 [lbs] (0.2361 %) ; CAPACITY = 1250.0 [lbs] 
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Fig. 4a Check load residuals without temperature correction. 
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Fig. 4b Check load residuals with temperature correction. 
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